Evolution Kernels of Twist-3 Light-Ray Operators 
in Polarized Deep Inelastic Scattering 

B. Geyer, D. Miiller, D. Robaschik 
Fakultat fur Physik und Geowissenschaft der Universitat Leipzig 
C ! Augustusplatz, D-04109 Leipzig, Germany 

i— s 

m 

7—1 ' Abstract 

> : 

The twist three contributions to the Q 2 -evolution of the spin-dependent 
structure function 52 (x) are considered in the non-local operator product ap- 
^sO \ proach. Starting from the perturbative expansion of the T-product of two elec- 

tromagnetic currents, we introduce the nonlocal light-cone expansion proved 
. by Anikin and Zavialov and determine the physical relevant set of light-ray 

operators of twist three. Using the equations of motion we show the equiv- 
O alence of these operators to the Shuryak-Vainshtein operators plus the mass 



operator, and we determine their evolution kernels using the light-cone gauge 
with the Leibbrandt-Mandelstam prescription. The result of Balitsky and 
'-^ ' Braun for the twist three evolution kernel (nonsinglet case) is confirmed. 



"*-! 1 Introduction 



According to the present status of quantum chromodynamics nonperturbative in- 
puts, being not calculable from first principles, are necessary for the interpretation 
and prediction of experimental results. Nevertheless, these inputs, for instance par- 
ton distribution and fragmentation functions, satisfy evolution equations whose in- 
tegral kernels are calculable perturbatively. For the twist two parton distribution 
functions all evolution kernels are now known up to next-to-leading order. 

However, also higher twist effects are experimentally accessible. Recently, in deep 
inelastic scattering the first moments of the polarized structure function g2{x) are 
measured |1] . In leading order of the momentum transfer Q 2 this structure function 
is determined by twist two as well as twist three contributions. In comparison with 
the twist two case the treatment of twist three is technically more subtle due to the 
appearance of a set of operators mixing under renormalization and constrained by 
relations between themselves. Up to now there exist already several papers which 
determine the local anomalous dimensions J| || |J as well as nonlocal operators and 
the evolution kernels for the distribution functions 0, |, [J. 



The two most complete calculations of the evolution kernels including the flavor 
singlet case, and 0, are based on complicated techniques; moreover, the compar- 
ison between both approaches is by no means trivial and both papers contain also 
some misprints. Therefore, in a first step, our aim is to confirm the known results 
using an independent approach and compare it with the local anomalous dimensions. 
In a next step we will study the solution of the evolution equations including the 
singlet case. Here we present the first results for the nonsinglet case. Technically 
we apply the nonlocal operator product expansion 0, || proved by Anikin and 
Zavialov in renormalized quantum field theory. To avoid the complications with 
ghost and gauge-variant operators which are present in covariant gauge we apply 



the light-cone gauge using the Leibbrand-Mandelstam prescription M, [10 



2 Perturbative Analysis 

The part of the hadronic tensor being relevant for polarized deep inelastic scattering 
is determined by the antisymmetric part (with respect to Lorentz indices) of the 
absorptive part of the virtual forward Compton amplitude. In the Bjorken region 
the leading terms in Q 2 of the Compton amplitude correspond to the leading light- 
cone singularities in the coordinate space. 

Our analysis starts with a perturbative investigation of the time ordered product 
of two electromagnetic currents. For simplicity, we will neglect flavor and color 
indices in the following. In leading order the antisymmetric part is determined by 
both twist two and twist three contributions. A perturbative expansion up to the 
first order in the coupling constant g provides additionally twist four contributions: 

{Tij(x)^(x)^y)YHy)} as = e^{id x p D c {x,y) i>(x) lal5 U{x,y)il;(y) (2.1) 
+ 9 -D c {x, y) ^(x) 7 .75 f o dt {2t - 1) U{x, z)F pX {z){x - y) x U(z, y)^{y) 

- ^D c (x, y) i){x)j a f dtU(x,z)F pX (z)(x - y) x U(z,y)ij(y)} + (x,fi <-» y,v), 
2 Jo 

where z = z(t) with z(l) = x and z(0) = y; F a/3 = l/2e a p fll/ F^ u is the dual field 
strength tensor and the path ordered phase factor 

U(x,y) = P exp j-z# J* dw^A^w)} (2.2) 

ensures gauge invariance. Instead of the free field propagator Sq(x) we used here 
the spinor propagator in an external gluon field (which is a straightforward gener- 
alization of the abelian case 01). One gets up to the first order in g: 

S c (x,0) = S c (x)U s (x,0) + 

9 -D\x) jf 1 dt {(At - 2)Yx u - Y tfl u } U s (x, tx)F^(tx)U s (tx, 0) + • • • , 

with 

U s {x,y) = Pexpj-i^ J dXA^wiX))^ - y^}, to (A) = x\ + y(l - A). 
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The nonlocal light-cone expansion is obtained by approximating the vector x by 
the light-like vector x defined as x = x + a(x, 7])r], where r\ denotes a fixed auxiliary 
vector (e.g. normalized by rj 2 = 1) and a the corresponding coefficient. In leading 
order we substitute px — > px, where p denotes a momentum, but the ^-singularities 
remain unchanged, x 2 — > x 2 . In general, the coefficient function will depend on 
two auxiliary variables Ki, whose range (according to the a- representation of the 
contributing Feynman diagrams) is restricted by < «j < 1. Here we introduce 
these variables quite trivially through integration over two 5-functions. In this more 
or less intuitive way we get the following light-cone expansion: 

{Tf(x)f(0)r = I' dm ['dK^ - l)<f(« 3 ) 

Jo Jo 

\id p D c (x) i)(K 1 x)^ a ^U s {K l x,K2x)i){K 2 x) (2.3) 

- Q -D c {x) V(^)7«x / 1 drf/ s (/ ?1 £,z)[7 5 (2r- l)F pX {z) 
I Jo 

- iF pX (z)}x X U s (z, K 2 x)tl)(K 2 x)} + («i ^ K 2 ), 

with z = [(k\ — k 2 )t + k 2 )]x. In renormalized quantum field theory the proof of this 
representation is more complicated (for a rigorous derivation see [0]). In the following 
we use the light-cone gauge xA = 0, so that U s (kix, k 2 x) = 1. Obviously, expression 
(p.3|) contains two light-ray operators from which we start our consideration: 

O p {ki,k 2 ) = i'(K 1 x)'J p 'J 5 ijj(K 2 x), (2.4) 

0[ P(T ](ki, k 2 ) = -^(ki£)7o. jf dr[j 5 (2r - l)F pX (z) - iF pX (z)]x x ij;(K 2 x) 

-(p<-<7). (2.5) 

The operator O p (k\,k 2 ) contains twist two as well as twist three contributions, 
whereas 0[ p<T ](/«i, k 2 ) contains contributions of twist 3 and higher (because of their 
coefficient functions the last operator is power suppressed). 



3 Choice of Light-Ray Operators 

As mentioned before, the operator ( |2.4j ) has mixed twist, so we look for a decompo- 



sition into its parts of definite twist. For the local light-cone operators there exists 
a well-known procedure: The tensor structure of operators with definite dimension 
decomposes into irreducible representations of the (formal) symmetry group 0(4). 
For light-ray operators, however, one has to take into account towers of such irre- 
ducible representations. Contraction with the light-cone vector x projects onto the 
leading twist two piece, 

O tw2 (fti, k 2 ) = x p O p {kx, k 2 ) = ^{kix) fi'y 5 ip(K 2 x). (3.1) 

The usual local twist two operators with spin n follow according to (d/dn = xd) 

O^ 2 = (ixd^O^O, «) lK=0 = ■ ■ ■ x^O™^, (3.2) 
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where O*™ 5 '.«„ is symmetrised and traceless. 

In order to construct the twist three part we represent the operator ( |2.4j ) in terms 
of local operators, extract the twist three part of them, and express the result again 
in terms of nonlocal operators. It turns out that 

0* w3 (ki,K 2 ) = — i(«2 — Kl) / duuO p (Kl, K\U + k 2 u), (3.3) 

Jo 



O p (k 1} k 2 ) = / duifj(Kix)i'y[pX cr ]D a (u } Kix, k 2 x)7 ip{[ K i u + k>2u]x), (3.4) 

«/ 

where we introduced the notation D p {u, k±x, k 2 x) = d p K2i + igA p ([Kiu + K2u\x) and 
u — 1 — u. Analogous to the twist two case our result can easily be checked by 
comparing it with the local operators of definite twist. 

All operators with the same or lower twist can be mixed under renormalization. 
Therefore, we also have to take into account the following twist 3 operators: 



Sp(Ki,T,K 2 ) = ig^(Kix) ft iF ap (rx) ± 7 5 F Q , p (rx) x a ^(K 2 x). (3.5) 



These operators are the building blocks of the operator fl2.5|) , which is contained in 
the light-cone expansion. For the local case, the importance of the operators (p.5|) 



has been first observed by Shuryak and Vainshtein in fl2fl . They result from an 
application of the equations of motion to the local twist three operators. 

The same can be shown directly and more easily also for the nonlocal operators. 
Applying the equations of motion to the operator O p and using the relation A u (ux) — 
A u {vx) = /" drx^Fpyijx) which is valid in the light-cone gauge we get 

6p(Kl, K 2 ) = ~(K2 - Ki) f du Utf° M (u) + nf EM (u) - 2M p (K lj K X U + K 2 U) 
— J 

+ U + S p (Ki, K\U + K 2 U, K 2 ) + U ~S p (Kx, KxU + K 2 U, K 2 )J, (3.6) 

where fi EOM is an equation of motion operator, f2 REM contains residual trace terms 
(being proportional to x p ) and operators containing an overall derivative, and M p 
is the mass dependent operator defined by 

— % 

M p (kx,k 2 ) = m tp(K 1 )a ap x a 'j 5 ( y xD)(K2x)i!(K2x), a Q/3 = ~[j a ,jp). (3.7) 

Corresponding to Eq. ( |3~3"1) also the operator 0* w3 is build up of the Shuryak- 
Vainshtein operators S p and the mass operator M p . Because of the following 
properties the additional operators f2 EOM and f2 REM can be neglected: If the op- 
erator VL is sandwiched between physical states their matrix element vanishes. 
The trace terms in the remaining operators Q REM can be safely neglected and the 
obtained overall derivative vanishes in the forward case. 
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4 Evolution Kernel of the Light-Ray Operators 



Before we present our result for the forward evolution kernel of the nonsinglet twist 
three light-ray operators we will discuss shortly the mixing properties of these op- 
erators. If we take into account all possible mixing partners up to trace terms, then 
we get the following set of operators 

{d p , %, M p , nf M }. (4.1) 

Also the equation of motion operator Q^ OM possesses an anomalous dimension and 
will be mixed with the other operators. However, from the general renormalization 
properties of gauge invariant operators it is to be expected that the counter term 
of Q^ OM is only given by the operator itself |14[|. In fact, this was explicitly shown 



in H| and recently in J3|. Thus, the anomalous dimension matrix is triangular. 
Forming physical matrix elements the equation of motion operators as well as their 
anomalous dimensions drop out completely from the renormalization group equation. 
Furthermore, because of the constraint (|3.6|) we can eliminate the operator O p and 
are left with the operators ^S p and M p only. 

There exists an additional property of the operators ^Sp which guarantees that 
+ S P and ~S P do not mix under renormalization. The contraction with e aa p p x a x*P , 
where x* is a second light-cone vector with xx* = 1, transforms the operators S p 
into ± ± S r <7 . = ± ± S' cr 7 5 (here the vertices of ^HS^ and ^Sa differ by a 7 5 matrix). Since 
the anomalous dimensions of ^5^ and ^HS^ coincide the operators ^Sp are effectively 
eigenstates of the duality transformation and thus do not mix under renormalization. 

Using FeynCalc and some own subroutines for doing the momentum integration 
in light-cone gauge we calculated the pole part of the contributing one-loop diagrams 
of the operators S p and M p in space time dimension n = 4 — 2e (the same set of 
diagrams as given in ||). Applying //j- = —^Q-§^ + • • • we get the result: 



[i 2 — + S p (k x ,k 2 ) = ^ f dy f dz 
dfi z 47T Jo Jo 

{(2C F - C A ) [y5(z)+S p (-K X y, k 2 - K X y) - 2z + S"{k x - k 2 (1 - z), -K 2 y) 

+K(y, z^Ml -y) + K 2 y, k 2 {1 - z) + k x z) 
+C A (2(1 - z) + L(y, z)Ys p (k 1 - k 2 z, n 2 y) + L(y, z) + S p (K iy , k 2 - k x z) 
+2C F (1 - yf 5(1 - z) MP(k 2 z - k iV )} (4.2) 

V 2 -T^~ SP ( K i> K 2) = ~r~ [ dy j dz 

dfi z 47T Jo Jo 

{(2C F - C A ) y5(zYS p (K 1 - n 2 y, -K 2 y) - 2z~S p (-K 1 y, k 2 - ^(l - z)) 

+K(y, z)~S p (ki (1 - y) + K 2 y, k 2 (1 - z) + k x z) 



(2(1 - z) + L(y, z)) S p (niy, k 2 - K X z) + L(y, z) S p (k x - k 2 z, n 2 y) 
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+2C F {1 - zf 5(1 - y) M p {k 2 z - Kl y)}, 
where 



(4.3) 



K(y,z) 
L{y,z) -■ 



1 + 6{z) + 8(y) — 

y z 

2 

6(l-y- z) + 5{z" ' 



1-2/ 



1-2/ 



and 



M p {k) 



AC, 



An 



dy[l + y- 



1-2/ 



M"(Ky). 



(4.4) 



To condense the notation we used ± S p (ki, k 2 ) = ± S' p (k 1 , 0, k 2 ), M p (k) = M p (0, k), 
and the standard plus-prescription fulfilling J* cfy [...]+ = and / dydz[..] + = 0, 
respectively; Cp = 4/3 and = 3 are the usual Casimir operators of SU C (3). 

Since charge conjugation transforms + S p (ki,k 2 ) into ~S p (k 2 ,ki) also the evolu- 
tion kernels in Eq. ( |4.2fl and Eq. ( [4.3[ ) are related to each other. It is easy to see 
that our kernels satisfy this symmetry condition. Up to a different definition of ± S P 
(and one small misprint) our result coincides with that of Balitsky and Braun || 
restricted to the forward case and m = 0. 

Finally, we want to add some technical remarks. 

Our operators are not completely traceless. Therefore, terms proportional to 
x p have been omitted. Furthermore, the auxiliary pole of the gluon propagator is 
regularized by the Leibbrandt-Mandelstam prescription 

1 kx* 



(4.5) 



with the light-like vector x 
proved in one loop order, e.g. 



kx (kx)(kx*) + ie 
. There are different advantages of this prescription 



consistency of tensor integral relations and the validity 
of power counting ||, The investigation of local anomalous dimension shows 
that independent operators as well as special nonlocal operators may appear |10 



However, their anomalous dimensions decouple from the anomalous dimensions of 
gauge invariant operators ]I4|]. Also in our calculation a x*-dependent operator 
appears and, as expected, does not contribute to the physical sector. 

Let us point out shortly, that in our calculation different prescriptions of the 
auxiliar pole in the gluon propagator do only affect terms which are proportional to 
two ^-functions. In the final result this effect is compensated by the wave function 
renormalization. As an example, we consider a triangular diagram of the operator 
± 5' p with two external fermion lines and one external gluon line 



d^k- 



i((pi+k)xK-i — (p2+k)xK2) 



((pi + k) 2 — m 2 )((p 2 + k) 2 — m 2 )k 2 
The exponential can be rewritten as 



kx 



i(pix)Ki-(p2x)K2 



kx 




(4.6) 



(4.7) 
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The first term in the curly bracket is analytic in kx so that the typical problem of 
light-cone gauge occure only in the second term which, however, contributes to a 
trivial ^-structure only. 
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